Ordering effect of Coulomb interaction in ballistic double-ring systems 
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We study a model of two concentric onedimensional rings with incommensurate areas A\ and A2, 
in a constant magnetic field. The two rings are coupled by a nonhomogeneous inter-ring tunneling 
amplitude, which makes the one-particle spectrum chaotic. For noninteracting particles the energy 
of the many-body ground state and the first excited state exhibit random fluctuations characterized 
by the Wigner-Dyson statistics. In contrast, we show that the electron-electron interaction orders 
the magnetic field dependence of these quantities, forcing them to become periodic functions, with 
period oc l/(Ai + A2). In such a strongly correlated system the only possible source of disorder 
comes from charge fluctuations, which can be controlled by a tunable inter-ring gate voltage. 
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An outstanding issue in the physics of mesoscopic sys- 
tems is the interplay between disorder and chaos, and 
electron-electron (e-e) interaction. Experimentally this 
problem is relevant for the anomalously large persistent 
currents in disordered metallic rings Q, magnetocon- 
ductance fluctuations through ballistic quantum dots Q , 
Coulomb peak height fluctuations || and peak spacings 
fluctuations Q in disordered quantum dots. Interaction 
is expected to play a crucial role in few-electron systems 
of low dimensionality; the way in which it affects quan- 
tum interference in phase coherent mesoscopic systems is 
presently under intense investigation. 

With this motivation in mind, we consider here a sys- 
tem of two concentric ballistic rings with incommensurate 
areas A\ and A2, in a perpendicular constant magnetic 
field B. The phase accumulated going around each ring 
is (p p = BA v /<pQ, p = 1,2, where 4>q = ch/e is the flux 
quantum. We will assume that the two rings are coupled 
by a site- dependent tunneling matrix element. A parti- 
cle will pick up random phases during its complicated 
motion. Such a breaking of the translational invariance 
can make the system nonintegrable and the one-particle 
spectrum chaotic, even in absence of e-e interaction. 

Our aim is to investigate how strong interactions affect 
the energy fluctuations of the many-body ground-state 
and the first excited states as a function of the applied 
magnetic field. We will show that the interaction orders 
the flux dependence of the lowest energy levels, which in 
the noninteracting case are random functions described 
by the Wigner-Dyson statistics. In contrast to the case 
of a clean single-channel ring || , in the two-ring geome- 
try interaction plays an essential role and modifies qual- 
itatively the period and the amplitude of the persistent 
currents. 

We choose a lattice model of two concentric one- 
dimensional rings with M sites per ring, described by 
the following second quantized Hamiltonian 



V (\ T i,p- v 3,P'\)'' 
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i,j,P,p' 



The operators c\ and c i create and destroy a spinless 
electron at site i of ring p. (i,j = l,...,J\f are site in- 
dexes and p, p' = I, 2 are ring indexes). The term propor- 
tional to t describes the intra-ring hopping, affected by 
the magnetic field. The site-dependent amplitude £12(2) 
represents inter-ring tunneling. The last term represents 
the translationally invariant e-e interaction. 

We first discuss the case of noninteracting particles. 
If the two rings are decoupled (£12^) = 0) the single- 
particle spectrum is given by e p (ki) = — 2icos[(fcj + 
2ttlp p /N)\ = -2tcos[(fe, + 2nH a p /N)], with k, = 
2-Ki/N, i = 0, ±1, . . . ± Af/2. Here we have introduced 
two independent, dimcnsionlcss parameters, H (mag- 
netic field) and a p (area of ring p). When H = 1 and 



= 1 the flux enclosed by ring p is 



Throughout 



n 



i,p i 



H.c. 



this paper we will consider two incommensurate areas 
di = [\/(5) — l]/2 ~ 0.618 and a 2 — 1.5. In that case 
the dependence of the single particle spectrum on H will 
be effectively random, with very complicated level cross- 
ings. Suppose now that the two rings are coupled by a 
tunneling amplitude t±2(i), small compared to the mean 
level spacing for each single ring, A = 8t/Af. To first or- 
der the tunneling will couple states of different rings. If 
ti2(i) is also site independent, the spectrum will hardly 
be modified and, in particular, level crossings will not 
disappear since they always involve levels of different mo- 
menta which are not coupled by a homogeneous matrix 
element. However, if t 12 (i) is not homogeneous, states 
of different momenta can be coupled, causing the lifting 
of the degeneracies. Such level repulsion is one of the 
essential characteristics of the Wigner-Dyson statistics. 
The single-particle energies can be made interact more 
strongly by increasing the magnitude of £12(2). We ex- 
pect that when the level splitting becomes of the order 
of A, a fully chaotic spectrum will develop. 
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FIG. 1. Probability distribution, P(s), of the one-particle 
energy-level spacings for a system of two rings, coupled by 
the tunneling amplitude £12(1) = —t/2 * cos 2 (37ri/27V)- Both 
rings have 11 sites; their area is ai — [a/(5) — l]/2 ~ 0.618 
and a,2 = 1.5 respectively. The probability density P(s) is 
obtained by performing a spectral average over the energies 
and an average over the flux. The dashed line represents the 
Wigner surmise plotted as a comparison. 

In order to show that the one-particle spectrum is 
chaotic we study its level statistics. Statistical proper- 
ties are obtained by performing a spectral average over 
many energy levels for fixed H and a subsequent aver- 
age over H . We consider the probability distribution 
of the nearest-neighbor level spacings, P(s), which we 
compute numerically by diagonalizing exactly a system 
of two rings, both having N = 11 sites, coupled by 
t 12 (i) = -t/2* cos 2 (37ri/2A/"). For each value of H the 
corresponding energy spectrum is first unfolded to ensure 
that the mean level spacing is unity |J . The distribution 
P(s) is then computed and finally averaged over many 
values of H . In Fig. 1 we plot P(s) averaged over 500 
values of H between and 5. The numerical curve is 
in good agreement with the Wigner surmise of the or- 
thogonal ensemble characterizing the Wigner-Dyson 
statistics. 

Now we might ask what are the consequences of such a 
spectrum on the properties of the noninteracting many- 
body system. The H dependence of an experimentally 
relevant quantity such as the total magnetic moment, 
m(H) — 8Eq (H)/dH, where Eq is the ground-state en- 
ergy for a N particle system, is just a very irregular func- 
tion of H , see Fig 2(a). However the chaotic nature of the 
one-particle spectrum reveals itself in the statistical fluc- 
tuations of the many-body ground-state and quasiparti- 
cle excitations as function of H and N. We are interested 
in two quantities: a) The distribution of the first excita- 
tion energy above the Fermi energy, SE = Ef — E^; b) 
The distribution of ground-state inverse compressibility: 



k = d 2 E<?/dN 
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Here E( 



N 



is the energy of the first excited state above the Fermi 
energy. It is obvious that for noninteracting particles at 
zero temperature 5E N and k n are the same and they are 
equal to £at+i — ejv, where en is the N th single-particle 
energy. Thus the fluctuations of 6E and re, averaged over 
the Fermi energy (namely over N) and H, scale with 



the A and their distribution is given by the same P(s) 
of the one-particle spectrum which, for our coupled-ring 
system, is given by the Wigner surmise. 

How is this complex scenario modified by strong e-e in- 
teractions? To answer this question we have resorted to 
exact numerical calculations on small systems, using the 
Lanczos algorithm to compute the fluctuations of many- 
body ground state and the first low-lying excitations as 
a function of Ep and H . We have considered systems 
with N = 11 sites per ring and N — 4,5,6,7,8 parti- 
cles. The number of electrons in each ring, N p , is not a 
good quantum number if the rings are coupled. The ex- 
pectation values (N p ) will however satisfy the constraint 
(N%) + (N 2 ) = N. We will first discuss the case of a short 
range interaction, coupling two nearest-neighbor sites 
(«,p) and (i + l,p) in the same ring with matrix element 
Vn, and opposite sites (i, 1) and (i, 2) in different rings 
with matrix element V12. The relevant energy parame- 
ter, indicating the relative strength of the interaction, is 
Vpp'/t. It is interesting to see how m{H) is modified. A 
plot of m{H) vs H for a system of N = 8 particles inter- 
acting with V\\ = 6.5 i, V12 = 6.3 1 is shown in Fig. 2a. 
The effect of the interaction is visibly spectacular. We 
can see that the aperiodic independence of m(H) of 
the noninteracting case is replaced by a perfectly peri- 
odic function, oscillating with period Hq — l/(a± + a 2 ), 
exactly as we would have for a single ring with area 
atot = (til +02)- The periodic m(H) develops smoothly 
with increasing interaction and eventually it becomes in- 
dependent of the interaction strength. The amplitude of 
m(H) in the strongly interacting regime is approximately 
one half of the amplitude of the noninteracting case. 
2 




FIG. 2. (a) Magnetic moment m(H) vs H, and (b), (c) 
four lowest many-body energies vs H for the two-ring sys- 
tem of Fig. 1. The total number of particles is N — 8. (a) 
The periodic solid line is the interacting case. The period is 
Ho ~ l/(ai + 0,2); the irregular dashed line is the noninter- 
acting case, (b) Noninteracting case, (c) Interacting case. 
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Moreover, not only is the ground-state energy a peri- 
odic function of H, but so are the first low-lying excited 
states, as we show in Fig. 2(c). In particular the low 
energy sector of the many-body spectrum is identical to 
the single-particle spectrum in one ring. Similar results 
are obtained for other values of even N. By computing 
the ground-state wavefunction we can also evaluate the 
particle density in each ring. The calculation shows that 
in the strongly interacting case (Ni) ss (N 2 ) ~ N/2 for 
any value of H . For odd N the interaction again forces 
the system to develop a periodic low-energy spectrum 
with the same period, l/(di + £12)- The main difference 
with respect to even N is that the H dependence of the 
spectrum is shifted by Hq/2. 

Strong intra-ring interaction is known to generate 
Wigner-crystal-like ground states in each single ring 
It has also been shown that inter-ring interaction can 
cause a nondissipative Coulomb drag |)| that affects qual- 
itatively the persistent currents in each ring |]Io| . Our 
exact results suggest that strong inter-ring interaction 
glues together the two stiff Wigner crystals created by 
the intra-ring interaction in each ring, resulting in a rigid 
structure that rotates like a single solid body under the 
effect of the magnetic field. Based on this picture, we can 
easily see why the oscillation period of the total magnetic 
moment of such a structure is l/(a\ +02). Indeed the 
Hamiltonian of the two frozen crystals glued together, 
each containing N/2 particles, is found to be 

a.2) 



H =2I 



d N . 



(2) 



where —i h-gg is the canonical momentum conjugated to 
6, the rotation angle of the solid around an axis passing 
through the center of the rings. I = m e ^{A\ + A%) 
is the moment of inertia. In zero flux, the eigenstates 
^n(0) = exp(in9) must satisfy the boundary condition 
\fr n (Q) = $ n (0 + 2tt/(JV/2)) due to the ordered crys- 
tal structure. This immediately implies the announced 
l/(ai + 0S2) periodicity as a function of H. 

We have also studied the case of an unscreened 
long range interaction coupling all the sites with ma- 
trix element V(\r, 



hP 



2^/7^7 



= V c R h 
7y cos(27r(i — j)/Af) 



— r 



3,P 



1/2 



The 



dimensionless interaction constant is a = V c R\/Tiv-f = 
(V c /t)[Af /47T sin(fcp)]. In this case we have to face first 
a new occurrence: in the strongly interacting regime 
(a w 5) where spectrum ordering should set in, the elec- 
tronic configuration corresponding to the minimum en- 
ergy is always the one with the inner ring empty and 
all the particles located in the outer ring. The corre- 
sponding low energy sector of the spectrum is obviously 
a periodic function of H with period l/a.2. To by-pass 
this unwanted situation we can add a term in the Hamil- 
tonian in the form of a potential difference between the 
rings, e.g. provided by a gate voltage: H g = £V p e p m t p. 
Typically a voltage difference V g = ei — 62 ~ —V c /2 en- 
sures an almost equal population in the two rings. In 



Fig. 3 we plot m(H) for N = 6 particles interacting 
with a Coulomb force of intensity V c = 10 1, for two val- 
ues of the voltage V g . Let us consider first zero inter- 
ring tunneling (solid line) (ll). We can see that, for 
V g = — 0.75 V c , m(H) is perfectly periodic with the ex- 
pected period l/(ai +02). On the other hand imperfec- 
tions and small glitches in the periodic pattern appear 
for V g = — 0.5 V c - In fact the calculation of the ring 
occupancy shows that while for V g — — 0.75 V c the elec- 
tronic configuration (Ni = 3,N 2 — 3) is stable for all 
values of H, when V g — 0.5 V c the two configurations 
(Ni = 3, N 2 = 3) and (N r = 2, N 2 = 4) alternate in the 
ground state. The aperiodic jumps in Fig. 3(a) corre- 
spond to a switch from one configuration to another as 
H varies. It is instructive to look at the H dependence 
of the lowest four energy levels, plotted in Fig. 4, and 
their corresponding ring density. Fig. 4(a) shows that 
when Vg = —0.5 V c these four levels are often two by two 
degenerate, the degeneracy being caused by the presence 
of the two electronic configurations [N\ — 3,N 2 = 3) 
and (Ni — 2, N 2 — 4) mentioned above. The gate volt- 
age is not yet strong enough to enforce a fixed electronic 
configuration in the two rings. As a result the levels 
are not completely periodic and they are sensitive to a 
inter-ring tunneling- see dashed line of Fig 3(a). How- 
ever if we increase V g up to 0.75 V c , the configuration 
(Ni = 3, N 2 = 3) is locked in the first three states. It 
is only starting from the fourth level that charge redis- 
tribution, with the appearance of the (N\ = 2, N 2 = 4) 
configuration at some values of H, is again possible. The 
important consequence is the ordering of the first three 
levels, which is robust against a nonhomogeneous tunnel- 
ing between the rings, as shown in Fig. 4(c) and Fig. 3(b). 
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FIG. 3. Magnetic moment vs H for the system of Fig. 1. 
There are N = 6 particles interacting with a long-range 
Coulomb potential of intensity V c = 10 1. The periodic solid 
line represents the case of £12(2) = 0. The oscillating pe- 
riod is l/(oi +0,2). The dashed line represents the case of 
ti2(i) = —t/2 cos 2 (27ri/jV). A voltage difference V g is applied 
between the two rings, a) V g = —0.5 V c . (b) V g = —0.75 V c . 
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J 12 *0 V g = -0.75V ' c 




FIG. 4. Energy vs H of the ground state and the first three 
excited states for the system of Fig. 3. (a) The two rings are 
decoupled \ti% = 0] and the voltage difference V g between the 
rings is -0.5 of the Coulomb interaction intensity, V c = Wt. 
(b) ti2 = and V 3 = -0.75 V c . (c) t 12 (i) = -t/2 cos 2 (2%i/Af) 
and V a = -0.7514. 



Notice that in the interacting case the fluctuations of k 
scale with the Coulomb energy. To obtain the probability 
distribution of 5E and k, for each value of H we first 
need to unfold the sequence of SE N and k n obtained 
upon varying N, to take into account the monotic change 
of Ajv- Once this is done, no fluctuations around the 
average are left and the distributions of 8E and n are 
just delta functions centered at s = 1 in rescaled units. 
Since this result is valid for any value of H, the average 
over H does not add anything. This is in sharp contrast 
to the case of noninteracting particles, where we saw that 
the two distributions were given by the Wigner surmise. 

In conclusion, we have shown that e-e interaction or- 
ders the flux dependence of the low energy spectrum of 
a two-ring system, suppressing the irregular behavior of 
the free particle case. The fluctuations of the ground- 
state and first excited state are no longer described by 
the Wigner-Dyson statistics. When strong e-e correla- 
tions are present, disorder can only be induced by charge 
redistribution controllable through an external gate. 

We thank I. V. Krive for useful discussions. Financial 
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< H < Hq/2 can be approximately parametrized by : 

E»(H) = ^N 2 +A N H\ (3) 

E?{H) = ^N 2 + A N (H - H) 2 , (4) 

for even N, while the expression for odd N are obtained 
shifting H by H /2. Here A N ~ A(N/Af) is a slowly 
increasing function of N and C is a constant. Thus for 
fixed and H we obtain, 

SE N (H) = A N (H 2 -2H H), (5) 
K N (H) =2%- A N (H 2 /2-2H H) (6) 
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